Search for: All records

This article presents an ultraweak discontinuous Petrov-Galerkin (DPG) formulation of the time-harmonic Maxwell equations for the vectorial envelope of the electromagnetic field in a weakly-guiding multi-mode fiber waveguide. This formulation is derived using an envelope ansatz for the vector-valued electric and magnetic field components, factoring out an oscillatory term of exp(-ikz) with a user-defined wavenumber k, where z is the longitudinal fiber axis and field propagation direction. The resulting formulation is a modified system of the time-harmonic Maxwell equations for the vectorial envelope of the propagating field. This envelope is less oscillatory in the z-direction than the original field, so that it can be more efficiently discretized and computed, enabling solutions to the vectorial DPG Maxwell system in fibers that are 1000x longer than previously possible. Different approaches for incorporating a perfectly matched layer for absorbing the outgoing wave modes at the fiber end are derived and compared numerically. The resulting formulation is used to solve a 3D Maxwell model of an ytterbium-doped active gain fiber amplifier, coupled with the heat equation for including thermal effects. The nonlinear model is then used to simulate thermally-induced transverse mode instability (TMI). The numerical experiments demonstrate that it is computationally feasible to perform simulations and analysis of real-length optical fiber laser amplifiers using discretizations of the full vectorial time-harmonic Maxwell equations. The approach promises a new high-fidelity methodology for analyzing TMI in high-power fiber laser systems and is extendable to including other nonlinearities. 

A prime objective of modeling optical fibers is capturing mode confinement losses correctly. This paper demonstrates that specific modeling choices, especially regarding the outer fiber cladding regions and the placement of the computational boundary, have significant impacts on the calculated mode losses. This sensitivity of the computed mode losses is especially high for microstructure fibers that do not guide light by total internal reflection. Our results illustrate that one can obtain disparate mode confinement loss profiles for the same optical fiber design simply by moving the boundary to a new material region. We conclude with new recommendations for how to better model these losses.